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In this work, we propose a method for calculating the free energy of anisotropic classical spin 
systems. We use a Hubbard-Stratonovich transformation to express the partition function of a 
generic bilinear super-exchange Hamiltonian in terms of a functional integral over classical time- 
independent fields. As an example, we consider an anisotropic spin-exchange Hamiltonian on the 
cubic lattice as is found for compounds with strongly correlated electrons in multiorbital bands and 
subject to strong spin-orbit interaction. We calculate the contribution of Gaussian spin fluctuations 
to the free energy. While the mean-field solution of ordered states for such systems usually has 
full rotational symmetry, we show here that the fluctuations lead to a pinning of the spontaneous 
magnetization along some preferred direction of the lattice. 


I. INTRODUCTION 

Recent research activities on transition metal oxides 
suggest that the interplay of the strong spin-orbit cou¬ 
pling (SOC), crystal field (CF) interactions, and elec¬ 
tron correlations may lead to compasslike anisotropic in¬ 
teractions between magnetic degrees of freedom.lil These 
anisotropic interactions have a generic form in 

which a depends on the direction of the particular link or 
bond and S denotes spin or pseudospin degrees of free¬ 
dom describing magnetic or orbital degrees of freedom. 

The models in which compasslike anisotropies are dom¬ 
inating, or also the pure compass models, have been 
known for a long time. These models appear naturally 
in strongly correlated electron systems as minimal mod¬ 
els to account for interactions between pseudospins de¬ 
scribing orbital degrees of freedom.l^H^ The compass-like 
anisotropies also arise as interactions between magnetic 
degrees of freedom in systems with strong SOC, which 
might be realized in 4d and 5d transition metal oxides.^ 
However, in these systems, due to the extended nature of 
4d and 5d orbitals, the compass interactions are always 
accompanied by the usual SU(2) symmetric Heisenberg- 
type exchange. These models are especially interesting 
because while the pure compasslike models are rare, the 
combined Heisenberg-compass models have been shown 
to be minimal models describing the magnetic properties 
of various materials. A review of the different realizations 
of compass models,f2Hll] their physical motivations, sym¬ 
metries, unconventional orderings and excitations may 
be found in the recent paper by Nussinov and van den 
Brink.E 

One of the common features induced by compasslike 
anisotropies is frustration, arising from a competition of 
interactions along different directions and leading to the 
macroscopic degeneracy of the classical ground state and 
in addition to rich quantum behavior. In many cases, the 
pure compass models do not show conventional magnetic 
ordering because the degeneracy of the classical ground 
state is connected to discrete sliding symmetries of the 
model.l^ Because these symmetries are intrinsic symme¬ 


tries of the model, they can not be lifted by the order-by 
disorder mechanisms. Instead, the direct consequence 
of the existence of these symmetries is that the natural 
order parameters for pure compass models are nematic, 
which are invariant under discrete sliding symmetries. 

The nematic order present in the compass model is 
fragile and is easily destroyed by the presence of the 
isotropic Heisenberg interaction which breaks some of 
the intrinsic symmetries of the model. In Heisenberg- 
compass models, some of the degeneracies become acci¬ 
dental. In these models, the true magnetic order might be 
selected by fluctuations via an order by disorder mecha¬ 
nism, removing accidental degeneracies and determining 
both the nature and the direction of the order parameter. 
Despite the simplicity of these models, the interplay of 
the Heisenberg and compass interaction leads to very rich 
phase diagrams even in the simplest case of the square 
lattice.l^ For classical systems this mechanism requires 
finite temperatures, where entropic contributions of fluc¬ 
tuations to the free energy become effective. 

In this work, we will be interested in studying the di¬ 
rectional ordering transitions in the Heisenberg-compass 
model on the cubic lattice.l^ From a historical perspec¬ 
tive, the three-dimensional 90°-compass model was the 
first model of this kind proposed by Kugel and KhomskiP 
in the context of the ordering of the t 2 g orbitals in tran¬ 
sition metal oxides with perovskite structure and then 
studied in more details by KhaliulliiJi^ in application to 
LaTiOa. The formal procedure which we will be using 
here is based on the derivation of the fluctuational part 
of the free energy by integrating out the Gaussian fluctu¬ 
ations, and determining which orientations of the vector 
order parameter correspond to the free energy minimum. 
To do so, we first express the partition function as a func¬ 
tional integral over classical fields. In this first paper, we 
consider classical spins at finite temperature. Our start¬ 
ing point in evaluating this exact representation of the 
partition function is the mean-field solution, which usu¬ 
ally does not reflect the anisotropic character of the inter¬ 
action referring to the crystal lattice axes. As a next step, 
we evaluate the contribution of Gaussian fluctuations to 
the free energy of the mean field ordered state. The latter 
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carries the information embodied in the anisotropic spin 
interaction and therefore allows to define preferred direc¬ 
tions of the spin order with respect to the lattice. We 
will not go beyond the simple evaluation of the contribu¬ 
tion of fluctuations, e.g., by incorporating the fluctuation 
contribution self-consistently. 

For simplicity, we choose the parameters of the model 
such that the ground state is ferromagnetic, i.e. we con¬ 
sider the Heisenberg interaction to be ferromagnetic and 
allow the compass interaction to be both ferromagnetic 
and antiferromagnetic. For any ferromagnetic and weak 
antiferromagnetic compass interactions, the minima of 
the fluctuational part of the free energy are attained if 
the spontaneous magnetization vector points along one 
of the cubic axes. 

This paper is organized as follows. In section II we 
introduce the functional integral representation of the 
partition function for the spin systems with interactions 
described by the most general bilinear form of the super¬ 
exchange Hamiltonian. The details of the method are 
outlined in the Appendix. In Sec. HI, we apply this 
framework to compute the angular dependence of the 
fluctuational part of the free energy for the ferromagnetic 
Heisenberg-compass model on the cubic lattice. Our re¬ 
sults are presented and discussed in Section IV. 


II. REPRESENTATION OF THE PARTITION 
FUNCTION 

We consider a system of identical classical spins S on a 
lattice, interacting in an anisotropic fashion as indicated 
in the introduction, defined by the Hamiltonian 


The partition function of the system is given by the 
integral over the Boltzmann weights of configurations 

Z = J[dS,] exp[-/3 5] - 1), (2) 

where /3 = l/fc^T is the inverse temperature, 5"“ are the 
components of the spin operator at site j. 

It is useful to represent the Hamiltonian in the basis 
of the eigenfunctions Xn;j,a of the spin exchange matrix, 
defined by 


E 


'^jj' Xn;j',ct' — f^nXn;j,a ■ 


For spins on a periodic lattice these eigenstates are la¬ 
beled by a wavevector q (inside the first Brillouin zone) 
and index v, characterizing three principle axes of the 
matrix J. Thus \n) = |q, and the normalized eigen¬ 
functions take the form 




iq R 


where N is the number of lattice sites, the u,y^a are or¬ 
thonormal real-valued eigenvectors, i.e., ,a = 

and Kq.j/ are the eigenvalues of the spin exchange 
interaction matrix. 

We now define the normal amplitudes of the spins as 


‘S'q.jy — ^ ' X<i,v,j,aSj 


and express the Hamiltonian as 


" = ■'S' ( 1 ) 

j,j' aa' 

where j,f label the lattice sites, a, a' = x,y,z label the 
three components of the spin and = 1. For the mod¬ 
els with compasslike anisotropic and Heisenberg isotropic 
interactions of spins, the interaction is diagonal in spin 
space, a = a'. The J“^“-matrix elements are different 
for the (j,j')-bonds along direction 7 with 'y = a and 
X a. However, since our consideration is also valid for 
the case when a ^ a', in the following, we will keep both 
indices. 

We will be interested in the long-range ordered phases 
of the system. The mean field approximation of the order 
parameter usually leads to a highly degenerate manifold 
of states, e.g., a ferromagnetic state with spontaneous 
magnetization pointing in any direction. This degener¬ 
acy is lifted by the anisotropic components of the spin 
interaction, but only at the level of the fluctuation con¬ 
tribution to the free energy (action) Sji. In the follow¬ 
ing, we outline a method allowing to calculate Sfi, which 
is based on the Hubbard-Stratonovich transformation of 
the partition function for spin systems described by the 
generic Hamiltonian 0 . We present details and discuss 
justifications for this method in the Appendix. 


H = ( 3 ) 

q.i^ 

where = S-q^,y. Commutation of classical spins al¬ 
lows us to employ a Hubbard-Stratonovich transforma¬ 
tion in terms of classical fields in order to represent 
the interaction operator as a Zeeman energy operator 
of spins in a spatially varying magnetic field. As a re¬ 
sult, one finds the following representation of the parti¬ 
tion function: 

Z = J [dip] (4) 

q,y 

where the integration volume element is given by 

iPdip*^„dipq^y 
27r|/tq,^ 

The contribution to the action in the case of classical 
spins is given by 

Sioci{<P*q,^,‘Pq,^}) = /3"^ y]ln[sinh(2^7)j)/2^(/jj], (5) 

3 
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where ip'j = , with (/?“ = 

Eq,,. and s(Kq,^) = 1 for < 0 

and s(Kq,,y) = i for Kq_^ > 0. The Hubbard-Stratonovich 
identity used to derive the above functional integral is 
different for eigenmodes with positive or negative 
eigenvalue Kq,iy, leading to the appearance of a complex- 
valued ipj. The details of evaluating Sioc 
can be found in the Appendix. 


III. APPLICATION TO THE CUBIC LATTICE 
A. Isotropic Heisenberg interaction 



(a) 


In order to demonstrate how to perform the evalua¬ 
tion of the above representation of the partition function, 
we consider first the isotropic ferromagnetic Heisenberg 
model with nearest neighbor interactions on the cubic 
lattice. In this case, the Hamiltonian Q reads 

^ = ( 6 ) 


where the lattice summation is over nearest neighbors 
(j,/)—bonds and J < 0. For the isotropic exchange 
interaction, the eigenvalues, Kq^^ = J^^cosq^, are in¬ 
dependent of V, Kq^i, = Kq, and hence are degenerate. 

A uniform ferromagnetic mean-field solution is found 
by solving the saddle point equation 


difMF difMF 


|Kq=o| ^{PA-IfY 


( 7 ) 


- /? ln[sinh(2,5(pMF)/2/3v?M^’] 


= 0 , 


where we used = '/NpMFSqfimo^^, pj = pmf, 

mo,i/ for the components of the unit vector along the mag¬ 
netization in the reference frame defined by the principal 
axes of the interaction matrix (which are the cubic axes 
in this case), and N is the number of lattice sites. The 
solution of Eq. 0 gives us a non-linear equation for the 
mean-field parameter: 

2|Kq=orV mf - 2coth(2/3(^MF) + -5-^^— = 0 (8) 

PPmf 

We solve this equation numerically at each temperature 
and get pmf{T). Linearizing Eq. ([^ near the tran¬ 
sition, we find the transition temperature = 

2|Kq=o|/3. We note in passing that a different length of 
the classical spin vector |S| = Sq may be simply scaled 
back to the unit length by changing the temperature as 
T' = SqT. Choosing Sq = 3/4 appropriate for quantum 
spin S = 1/2, we find the renormalized transition tem¬ 
perature T/ = |Kq=o|/2 , which agrees with the quantum 
mean-field transition temperature. 

The fluctuation contribution is obtained by expanding 
the action in the fluctuation field 6pq^,j = pq^i, — Pq Y 


FIG. 1: (Colors online) The magnitude of the action iSfl(0, (j>) 
dehned by Eq. ( |18[ ) is plotted on the surface of the unit sphere. 
The preferred directions of the magnetization, corresponding 
to the minima of the free energy, are shown by deep blue 
color. The energy scale is shown in units of J. J = — 1 and 
K = 0.75: the preferred directions of the magnetization are 
along the cubic axes. 


about the mean field solution to the lowest order: 


S = iSo -t- (9) 

Sq = N\Kq=o\~'^pY[p - NT\n[sivL\l{2j3pMF)/‘iYpMF\■ 

¥ov Gaussian fluctuations, the fluctuation part of the 
free energy, or equivalently the action, iSfi, is a bilinear 
function of 5pq^^. It is given by 


^ Aq^i,^>5p*q,,5pqy, ( 10 ) 


where we defined matrix elements of Aq ,ji,i describing 
the weight of the Gaussian fluctuations of wavevector q 
and polarization v as 


A, 


1-1 


5vy - F 


3 L' 


Pc{5^y - mo^^moy) + 


3l3rmo^„moy]s{Kq^„)s{Kqy) 


( 11 ) 


Here, for shortness we introduced r = l/{2j3pMFY ~ 

1 /sinh^(2,5<pMF)- 

In the limit of small q, it is instructive to separate 
the fluctuations into longitudinal (along mo) and trans¬ 
verse (perpendicular to mo) components, = uiQ-Spq 
and dt/Sq = 2 respectively. We defined 

“ (moXz)/| sin0| and 
m 2 = mi X mo, where cosO — mo-z. Despite the com¬ 
plex nature of fluctuational fields, their separation into 
transverse and longitudinal modes is possible in the limit 
of small q, because the interaction eigenvalues Kq i, < 0 
and thus, s{Kq^i,) = 1 in this region of the BZ for any 
polarization component f. Then, the longitudinal fluc¬ 
tuations contribute to the free energy as 

= E ("ioA¥>q)(mo-V_q). ( 12 ) 

Q 
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The transverse fluctuations are gapless in agreement with 
Goldstone’s theorem: 




q,Kq <0 




(13) 


since limq^o [|Kq| ^ - |/3c] = 0. 


B. Fluctuations due to anisotropic compass 
interactions 

Next, in addition to the isotropic Heisenberg term, 
let us take into consideration an anisotropic compass 
interaction, K. The constraint that the ferromagnetic 
mean field solution remains stable is satisfied for all neg¬ 
ative (ferromagnetic) values of K and for positive values 

ki¬ 
ln the presence of the anisotropic compass interaction, 
the model Q reads 

(«) 


breaks rotational invariance, which results in a selection 
of preferred directions of the order parameter, which min¬ 
imize the free energy. 

The 3 X 3-matrix may be diagonalized and 

has eigenvalues and eigenvectors v.y q, 7 = 0 , 1 , 2 . 
This allows us to express = 

where Sifq^-f = v^^q ■ The in¬ 
tegration over the fluctuation amplitudes may now be 
performed and gives 

Sfi = l3-^lY.^n\Xo,qXi.,qX2,ql (18) 

q 


where we chose s{Kq.,y) = ±* for Kq.i, > 0 , following the 
procedure described at the end of the Appendix. Alter¬ 
natively, we may use that |Ao,qAi,qA 2 ,q| = | det{Aq_iqy/}|, 
saving the trouble of having to determine the eigenstates 
of . 

Let us now derive the explicit expression for the fluc¬ 
tuation contribution for an arbitrary orientation of mo = 
(sin0cos(/), sin0sin(^,cos0). Inserting this into the defi¬ 
nition of Aq given by Eq.( 11), we find its elements to 
be 


where the exchange interaction is given by 

Jjj, = -Sji j+r[J + K6a^\r\] (15) 

The index r = zLx,zLy,±z labels nearest neighbor sites, 
where |t| = x,y,z specifies a direction in spin space (x 
for bonds along the x-direction, etc.). The eigenvalues 
of the operator J““ defined in the previous section are 
given by 


^ A"(5cK^iy) cos (^Q.. (15) 

Oi 

The eigenvectors are again along the three cubic axes, 
such that the components are This time the 

three eigenvalues for given q are not degenerate (except 
in the limit q —>■ 0 ) and the fluctuation contribution to 
the free energy will therefore depend on the orientation 
of the spontaneous magnetization. We may again use the 
representation of the partition function Z as a functional 
integral over the Fourier components fq^i, of the auxil- 
liary field. 

Provided J < 0 and K < |J|, the mean-field solu¬ 
tion (fiMF is given as before by solving the transcendental 
equation Q numerically. The fluctuation contribution to 
the free energy is obtained by expanding the action in the 
fluctuation field about the mean field solution to lowest 
order. We get 

Z = C exp(-/35o) J [d6(p] exp(-, 55 fi{ 5 v 5 q,,.}), (17) 


where the fluctuation part of the action is given by 
Eqs.(lO) and (11) In the following, we show that by 
comparison to the isotropic model, Eq. 0 manifestly 


■^q,01 — 3 ) (3/3r 

^< 7,10 = ^ q ,01 

■^q,02 — ^ )'® (^< 7,2 ) (3/?r Pc)C(pC0S0 

^q,20 = ^q,02 

Ag,ll = \Kq,y\~^ - ls{Kq.y)s{Kq^y){l3c{l - 

■^q ,12 3 )'^(^< 7 , 2 ) (3/5^ Pc)s4)^6^0 

Aq,21 = Aq^l 2 

Aq^22 |s:q, 2 | ^^{,^q,z')^(,^q,z^(^Pc^0 T 8 ^rC^), 

(19) 

where, to shorten notations, we denote sind(^) = 
S 6 )( 0 ) and cosd(^) = cg( 0 ). The interactions are de¬ 
fined as = 1/[(J -I-AT) cosga; -I- Tcos^y-|-Jcosg 2 ], 
Kq y = 1/ [(J -I- AT) cos Qy + J cos Qx + J cos Qz] and 

Kql, = 1/ [(J -I- AT) cos Qz + J cos Qx + J COS Qy]. We see 
that the matrix Aq^^^i has a rather complex structure as 
a function of q and angles 6 and (j). This gives rise to a 
complex behavior of the eigenvalues Ao,q, Ai^q and A 2 ,q. 

IV. RESULTS AND DISCUSSIONS 

We now present the results obtained for Sfi{9,(j)) by 
performing numerical integration in Eq. The angu¬ 
lar dependence of Sii{9,(f>) is presented in Figs.l, where 
the magnitude of as a function of orientation 

of the spontaneous magnetization is shown as a color- 
coded plot on the unit sphere. The calculations in Figs.l 
are performed at temperature /3 = /3c + 1 and assuming 
J = —1. We see that Sfi{9,(j)) has a non-trivial depen¬ 
dence on the direction of the order parameter defined 
by angles 9 and (j). This peculiar angular dependence of 
iSfl( 0 , (/)) is inherited from non-trivial angular dependen¬ 
cies of Ao,q, Ai^q and A 2 ,q. 
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In Fig.l, we present the profile of Sfi{6,(f)) computed 
for K = 0.75. We can see that Sfi{6,(j>) is minimized 
when the magnetization is directed along one of the cu¬ 
bic axes. We note that the cubic directions are also se¬ 
lected for other values of the compass interactions, both 
antiferromagnetic as well as ferromagnetic, where the fer¬ 
romagnetic state is the mean field solution {K < | J|). 


V. CONCLUSION 


The magnetic properties of heavy transition metal ox¬ 
ides such as iridates and others are emerging as a new 
fascinating field offering opportunities to realize strongly 
frustrated quantum spin systems in the laboratory. In 
these systems, the combination of multiband electronic 
structure and strong Coulomb and Hund’s couplings with 
strong spin-orbit interaction can give rise to extremely 
anisotropic spin exchange interactions of the compass 
type. Mean field solutions of these models are often un¬ 
touched by the anisotropies of the model and show the 
full isotropy of pure Heisenberg models, in contrast with 
experimental observations. In this paper, we addressed 
the question how the system selects special preferred di¬ 
rections of the mean field order parameter vector. We 
restricted ourselves to the case of a ferromagnetic order 
parameter, but an analogous question exists for antifer¬ 
romagnetic or more complicated ordered structures. We 
find that the high degeneracy of the ferromagnetic mean- 
field solution is lifted by the free energy contribution 
from thermal fluctuations. We calculated the fluctuation 
contribution for a Heisenberg-compass model of classical 
spins on a three dimensional cubic lattice with nearest 
neighbor interactions - an isotropic Heisenberg coupling 
J < 0 (which we take as the energy unit), and a compass 
coupling K. The ferromagnetic state is found UK < | J|. 
Rather than exploring the full phase diagram, we focused 
on one typical temperature T = Tc/{1 + Tc), where 
is the mean-field transition temperature. For values of 
K < 1, the system is found to choose preferred direc¬ 
tions of the spontaneous magnetization along one of the 
cubic axes. In the temperature regime considered here, 
we expect the classical approximation to be valid. A 
generalization to quantum spin systems of the approach 
presented here is in preparation. 
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Appendix A: Hubbard-Stratonovich transformation 
of the partition function for spin systems 


1. General formulation 


The Hubbard-Stratonovich (H-S) transformation is 
based on the mathematical identitiy 


exp[— 



where we defined 


dyexp 



2s{a)xy) , 

(Al) 


s 



1 , 

G 


if a < 0 
if a > 0. 


(A2) 


For a > 0 we may as well use s{a) = —i. We will 
later make use of this ambiguity when we evaluate the 
y-integrals approximately, which may lead to imaginary- 
valued contributions. 

In the above H-S-transformation, x may be a number 
or an operator. In the case it is an operator, we use the 
eigenfunctions \n) of a; defined by 


a;|n) = a;„|n) 


to prove that 


exp —axP' 



\n) = exp 
dy exp - 

dy exp - 


[-axl] \n) 

-f 2s{a)xny) 

|a| 

1 ^ -f 2s{a)xy) 


\n) 

\n). 


(A3) 


This identity also works for complex (non-Hermitian) 
X and y: 


exp[—aai^a;] 


i 

27r|a| 


dy* dy exp 



-I- s(a)(a;^y -I- H.c.) 


We now turn to the case of the partition function of 
a spin system with generic interaction Hamiltonian (1). 
In order to use the mathematical identities we need to 
represent the Hamiltonian (1) in terms of normal coordi¬ 
nates. To this end we define the normalized eigenstates 
of the exchange interaction operator 


'y ' jjj' Xn-,j',a' — KnXn;j,aj 


j',a' 

in terms of which we have 


J??' = Y 


^nXn-J.aXn-J',a’ j 


(A4) 

(A5) 


where Xn-j'.a' form a complete and orthonormal set of 
eigenfunctions and thus obey 

y ', Xn-.i,aXn’;j,a = dn,n'i (-^ 6 ) 

y ', Xn;i,aXn-,j',a' = djjrda^a'- 
n 
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For spins on a periodic lattice, the eigenstates \n) = 
|q, v) are labeled by wavevector q and spin component 
V, and the eigenfunctions take the form 

(A7) 

where are normalized real valued eigenvectors, i.e. 
Sa the eigenvalues of the spin 

exchange operator. We now define the normal amplitudes 
of the spin operators as 

(A8) 

j,a 

and express the Hamiltonian (1) as 

H = (A9) 


where S'* ,, = S-q^^. 

We seek to apply the above mathematical identities 
(A1)-(A3) to each normal component separately. This 
requires the normal components of the spin operators 
to commute with each other, which is certainly true for 
the classical spins. Then using the Hubbard-Stratonovich 
transformation one may express the Boltzmann weight 
operator of each normal mode in terms of normal field 
amplitudes (pq^i, as 

6xp[—/3Kq^,ySq j,Sq^l/] = j" j j ^V^q,!/'^^q,l' 

(AlO) 

SXp [ /3{|Kq,j/| ^q,u^<i,l^ “f ^i^(l,L'){Sq i^ipq ,j + Ff.C.)}] 

The complete Boltzmann weight operator may be ex¬ 
pressed, again using the commutability of the normal 
mode operators, as 

exp[-^ ^ Kq,,.S* ,,Sq,i.] = Adi/j] (All) 

q,i^ 

Gxp[ /3 'y A (‘^q,l^‘dq,l^ T ^-C')}]) 

q.s^ 


where <fq^ = ■ The integration volume element is 

given by 


[dip] = Tlqi, 


27r|Kq,i,| 


Next, we find that the partition function of an interacting 
classical spin system on an infinite periodic lattice may 
be expressed as 


Z = 


= C y [dp] exp 


(A12) 

~P 'y ^ l^q,i^l ^q,i/^q.,v ~ Sloc{,{Pq,v}) 


q.i^ 


where C is a constant. The contribution iSioc({</5q,i/}) to 
the action is given by 

‘5ioc({‘/5q,!/}) = ^ Ill ^0 (A13) 


and Wj is computed by taking into account the constraint 
of the unit length of classical spins, S| = 1, and integrat¬ 
ing over all directions of spin at each lattice site: 


W,= 


dSj dflj 


2tt 


exp 


dfli 


exp 


J 47 r 

sinh 2/3]pj] 


oc 

2I3Y.VP: 


m 


2P\pj\ 

This gives 

c’ 1 v^,„rSinh2/3|(/jj-| 


2p\p, 


1 ) 

(A14) 


(A15) 


Here we defined the complex-valued three-component 
field p'j at each lattice site j as 


(A16) 

q,i^ 

q,i' 

= ^R,j + i^py 

Observing that Kq^^, = K-q.y, we get 

^R,j = R-e{(/3“} = 51 ^^,l'X*q,u- 3 ,a (A17) 

= 51 ^^,’^X*q,u-j,a- 


The field amplitude is determined by 


‘fj = \J{^R,j+i^i,iY^ (A 18 ) 

where = {p\p p%j, Pr^^ and = 

(^pIj, pi^, pi^y 


We now derive the contribution of Gaussian fluctua¬ 
tions to the free energy for the ferromagnetic mean field 
state which we denote as Pmf- To this end, we expand 
‘5ioc({‘/?q,i/}) ( A15 I in terms of the fluctuation amplitudes 
and separate the mean-held and huctuational contribu¬ 
tions. First, we expand the Held amplitude pj to bilinear 
order in the Huctuation amplitudes: 


Pj = Pmf + Spj, (A19) 

~ {‘Pmf ■ {^‘Prj 

^Pmf 
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Using Eq. (A17), we now obtain the expressions for 5ipj 
and Sif'^ in terms of and keeping quadratic 

(Gaussian) terms only: 


2ipMF 




(A20) 




Next, we expand Eq. (A15) step by step as 

sinh2/?|(^j| = sinh(2/3((/3MF + ^<Pj)) 

= sinh(2/3</5MF)[l + 2{j55ipj)‘^] + cosh{2 j3 ip mf)2 (5 5 

and further 

In [s\nh{2(3(f j)/2(3tpj] 

= In [sinh(2/3((/3MF + - In [2/3 ((/?mf + 5p}j)] 

= In [sinh(2/3(^MF)/(2/3:^MF)] 

Sipj 


+ [2PipMF COth{2/3ipMF) - 1 ] 


1 

+- 


{2i3(pmf)‘^ 


+ 1 


‘Pmf 

^5ipj_Y 


2 [ sinh^(2/3(^MF) J ‘Pmf 

The fluctuation part of the local part of the free energy 
is then given by 

-/3"^<5^1n[sinh(2^(Pj)/2/3v3j] = 


= -fPc^mf Y^ ^^3 


(A21) 


1 M (2f3ipMFY 


2f3p>MF sinh {2(3ipmf) 




where we have used that 2 j3ip mf coih.{2 [3i^mf) ~ 1 = 
\j3cl3ip\jp ■ Substituting the expressions for 6p)j, Sip'^ and 


defining r = \/{2l3ipMFY — l/sinh^(2/3(^MF), we get the 
fluctuation contribution to the free energy 

Ac^p^'5^p*^„5lpc^y, (A22) 


where we defined matrices describing the weight of 

Gaussian fluctuations of wavevector q and polarization u 
as 


Aq.i/iy' — |Fq.iy| (A23) 

2 

-^[Pc{5yy - mopuioy) + 3Prmopmoy]s{Ktip)s{Kqy) 

The fluctuation matrix Aq^i,' will in general be non- 
Hermitian, and its eigenvalues will be complex. We now 
use that Aq^^i is an even function of q and divide q-space 
into Qx > 0 (M>) and Qx < 0 {M<). Note that the choice 
of Qx for dividing the BZ in half is arbitrary, and we could 
also do it with a help of Qy and Qz- For modes <Pqp with 
q we choose s(Kq,i/) = +i, whereas for modes with 

q GM< we choose s{Kqp) = —i, where Kq,,^ > 0 in both 
cases. Then we have A_q_i,i,/ = A* and as a result of 
the functional integration we will get 


Z = Zmf j [d5<p\ exp 




= Zmf exp 


= Zmf exp 


~o Y ln(det(Aq,^,,Odet(A* 


qeM> 


'y ) In I det(Aq_i,[^/) 


(A24) 


where 

Zmf = exp [-/SA^lFq^orVl/F] [sinh(2/3(/?MF)/(2/3v5MF)]^ ■ 
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